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TITLE OF THE INVENTION 



DEFECT INSPECTION APPARATUS AND METHOD 

CROSS-REFERENCE TO A RELATED APPLICATION 
[0001] This application is based upon and claims the benefit of priority from the prior 
Japanese Patent Application No. P2000-297430, filed on September 28, 2000, the entire 
contents of which are incorporated herein by reference. 

BACKGROUND OF THE INVENTION 
FIELD OF THE INVENTION 

[0002] The present invention relates in general to a defect inspection apparatus and method 
for detecting defects in a substrate having more than one pattern. Examples of substrates 
having such patterns include, but are not limited to, photo-masks or reticles, liquid crystal 
displays (LCD) panels, and printed circuit boards (PCB). 

DISCUSSION OF THE BACKGROUND 

[0003] Large-scale integrated circuit (LSI) patterns for semiconductor devices are transferred 
onto wafers from photo-masks. The photo-masks are fabricated so that mask writers such as 
an electron beam exposure apparatus or a laser-beam writer apparatus can write a plurality of 
desired patterns. If defects such as pin-dots or pinholes are located on the photo-masks, then 
they can have a significant influence on the resultant circuit patterns that are transferred onto 
the wafers. This would cause critical problems, including, but not limited to, drastic 
reductions in performance and yield. 

[0004] To avoid these problems, it is essential to perform defect inspection and repairing, if 
necessary, after the photo-masks have been written by mask writers. To this end, the 
employment of a cutting-edge defect inspection apparatus is required, which detects such 
detrimental defects in the photo-mask patterns with 100 percent accuracy. For the generation 
of 130 nm nodes by photolithography, defects as small as 100 nm need to be detected. 
Additionally, in recent years, in order to further minimize the critical feature size of circuit 
patterns on wafers, advanced photo-masks employing resolution enhancement techniques 
have been adopted, such as phase-shift masks and optical proximity effect correction (OPC) 



masks. 

[0005] Typical examples of phase-shift masks for resolution enhancement are the attenuated 
phase-shift mask that uses as its "absorber" a halftone membrane such as molybdenum 
silicide (MoSi) or the like instead of chromium membrane, as in a conventional binary mask, 
5 and the alternating phase-shift mask, the so-called Levenson mask, which is fabricated with 
quartz glass engraved on a binary mask. Either one of these masks is such that both optical 
amplitude and phase are controlled precisely, which is a technique for improving 
photolithography margins, including the exposure dose and focal depth. 
[0006] On the other hand, the OPC masks are for suitably setting critical dimension at a 

10 desired value while at the same time precluding non-uniformity of such critical dimensions, 
which otherwise can vary in accordance with the pattern density of the circuit pattern to be 
transferred. This is the technique used for correcting a critical dimension that is variable in 
value due to exposure by adding an auxiliary pattern to, or subtracting an auxiliary pattern 
from, the original on-the-mask pattern size, to thereby correct non-uniformity of the resultant 

1 5 pattern sizes. 

[0007] A prior known approach to inspecting phase-shift masks and optical-proximity-effect 
correction masks is to employ the following methodology. Firstly, uniformly illuminate a 
photo-mask with light while simultaneously scanning the pattern by use of an X-Y stage. A 
charge-coupled device (CCD) or the like is used to photo -electrically scan such light- 
20 irradiated patterns. With this photo-mask image being used as scanned data, apply image 

processing thereto. Having prepared in advance, or generating simultaneously, reference data 
for use as a standard image, compare the reference data and the scanned data to identify the 
differing portions as defects. 

[0008] Currently available reference data generation methods typically include a die-to-die 
25 comparison scheme and a die-to-database comparison scheme. The die-to-die comparison 

scheme is such that the above-noted method is employed to scan the optical image of a photo- 
mask with the same pattern being formed thereon for test use, which is then subjected to 
image processing. Resultant image-processed data is adaptable for use as the reference data 
required. The die-to-database comparison scheme is such that the expansion data expanded 
30 from the design data of a circuit pattern is used as the reference data. 

[0009] The die-to-die comparison scheme assumes that the test-use photo-mask is virtually 
free from any defects in the same portion of the dies that are compared to each other; 
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unfortunately, this assumption is not justified: a few of defects can still be present. In cases 
where defects happen to exist commonly in a region to be compared to a photo-mask of the 
reference, such defects can possibly be overlooked. 

[0010] Alternatively the die-to-database scheme offers the capability to obtain ideal 
reference data due to the fact that the design data per se is subject to simulation for use as the 
reference data. However, there is a tradeoff between higher sensitivity and computation 
complexity concerning the fidelity of reference data to the scanned data. To improve the 
fidelity, simulation methodology of high-level precision must be used to obtain the intended 
reference data from the design data, which would result in an unwanted increase in the 
amount of calculation. Additionally, with prior art simulation methods, results fail to match 
real images in the case of phase- shift masks, due to the absence of any careful consideration 
about interferences due to possible optical phase differences. With the optical-proximity- 
effect correction masks, minimizing the size of an auxiliary pattern can rapidly increase to the 
extent that it becomes as small as one fourth of minimum feature size, a widely accepted 
trend in micro-fabrication using photolithography technologies. It has also become difficult 
to simulate such ultra-fine or "micro" patterns. 

SUMMARY OF THE INVENTION 
[0011] It is therefore an object of the present invention to provide a new and improved defect 
inspection apparatus and a defect inspection method capable of creating reference data 
through simulation of an accurate optical image from design data without suffering from any 
risks of erroneous recognition. 

[0012] In accordance with a first aspect of the invention, there is provided a defect inspection 
apparatus for inspecting defects by optically scanning a pattern on a test substrate having the 
pattern, converting it into scanned image data for use as electrical image information, and 
then comparing said scanned data with reference data indicative of an optical image 
obtainable from design data of said test substrate, including: 

a first processor unit creating binary or multi- value expansion data by expanding said 
design data and obtaining one of complex transmission distribution data and complex 
reflection distribution data of said test substrate; 

a second processor unit calculating said reference data by using one of said complex 
transmission distribution data and said complex reflection distribution data passed through 



multiple complex finite impulse response (FIR) filters; and 

a comparator unit comparing said scanned image data with said reference data. 
[0013] In accordance with a second aspect of the invention, there is provided a defect 
inspection method for inspecting defects by optically reading a pattern on a test substrate 
having the pattern, converting it into scanned image data for use as electrical image 
information, and then comparing said scanned image data with reference data indicative of an 
optical image obtainable from design data of said test substrate, including: 

creating binary or multi-value expansion data by expanding said design data; 

obtaining one of complex transmission distribution data and complex reflection 
distribution data of said test substrate; 

calculating said reference data by using one of said complex transmission distribution 
data and said complex reflection distribution data passed through complex coefficient 
multiple FIR filters; and 

comparing said scanned image data with said reference data. 

BRIEF DESCRIPTION OF THE DRAWINGS 
[0014] A more complete appreciation of the present invention and many of the attendant 
advantages thereof will be readily obtained as the same becomes better understood by 
reference to the following detailed description when considered in connection with the 
accompanying drawings, wherein: 

Figure 1 is a diagram schematically depicting the overall configuration of the defect 
inspection apparatus in accordance with one embodiment; 

Figure 2 is a block diagram to explain the data processing method implemented in the 
defect inspection apparatus in accordance with the embodiment; 

Figure 3 is a block diagram showing the reference data calculation method in the 
defect inspection apparatus in accordance with the embodiment; and 

Figures 4A and 4C depict optical images of bright field, and Figures 413 and 4D are 
optical images of dark, wherein Figures 4A-4B are optical images due to reference data 
obtained by filter calculation in accordance with the embodiment, whereas Figures 4C-4D are 
optical images as calculated by applying a strict partial coherent model to the same pattern. 



DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 
[0015] One embodiment will now be set forth in detail with reference to the accompanying 
drawings below. 

[0016] Referring to Fig. 1, Fig. 1 schematically shows an overall configuration of the defect 
inspection apparatus in accordance with the embodiment. 

[0017] A test substrate (also called a sample or specimen in some cases) 30, having a pattern, 
is stably situated by vacuum adsorption on an X-Y stage 3 1 of the defect inspection 
apparatus, shown herein. Typical examples of the test substrate 30 are substrates with a 
multi-layer pattern such as a photo-mask, a reticle, a liquid crystal display (LCD) panel or 
equivalents thereto. Assume here that the test substrate is a photo-mask. 
[0018] Light emitted from a light source 32, such as a high-pressure mercury lamp or the 
like, is reflected by a reflex mirror 34 to fall onto the test substrate 30. A condenser lens 33 
uniformly illuminates the light irradiated from the light source 32 onto the test substrate 30. 
An objective lens 35 causes light that has passed through the pattern on test substrate 30 to be 
focused onto the image sensitive surface of a charge-coupled device (CCD) sensor 36. 
[0019] At this time, the optical resolution is greatly improved, with a decrease in wavelength 
of the light source 32 and also with an increase in numeral aperture (NA) of the objective lens 
35. 

[0020] An optical image of the test substrate 30 is focused on the image sensitive surface of 
the CCD sensor 36. This optical image is converted by CCD sensor 36 into corresponding 
electrical image data. This image data is then digitized by an analog-to-digital (AD) data 
converter 37 into a digital scanned image data, which will be input to a comparator device to 
be compared with reference data, to be next explained below. 

[0021] Figure 2 is a block diagram that explains the data processing method used during 
defect inspection. 

[0022] The scanned image data, as output from the AD converter 37 of Fig. 1, is input to a 
scanned image data to the input unit 5, shown in Fig. 2, for temporal storage therein. At this 
time, design data is simultaneously input to the unit 6 to obtain complex amplitude 
transmission distribution data (or alternatively complex amplitude reflection distribution 
data). At this stage, design data of each layer of the multi-layer pattern is used to create 
binary or multi-level value data to calculate complex amplitude transmission distribution data 
(or complex amplitude reflection distribution data) of the test substrate 30. 



[0023] The processor unit obtains the complex transmission distribution data (or complex 
amplitude reflection distribution data) operates to generate a set of binary or multi- valued 
data through expanding each design data of the multi-layer patterns and then provide the 
complex transmission distribution data (alternatively reflection distribution data) that is 
obtained through multiplication of this expanded data by a complex amplitude transmission 
(or complex amplitude reflection) and summing the two multiplication results. 
[0024] The complex amplitude transmission distribution data (or complex amplitude 
reflection distribution data) thus calculated through this processing is input to a reference data 
calculation unit 7, which is for calculating the reference data to be compared. Preferably this 
reference data calculation unit is arranged to perform processing by letting the complex 
transmission distribution data (or complex reflection distribution data) pass through multiple 
complex FIR filters. 

[0025] This reference data, output from the reference data calculation unit 7, and the scanned 
image data output from the scanned image data input unit 1, are compared to each other at the 
comparator unit 8 to detect any possible defects. 

[0026] An explanation will next be given of processing of unit 6, to obtain the complex 
amplitude transmission distribution data, and unit 7, to calculate the reference data. As in the 
case of the complex amplitude reflection distribution data, similar results are obtainable by 
setting a physical property value at the complex amplitude reflection coefficient in the 
following explanation; thus, a repetitive explanation thereof will be eliminated herein. 
[0027] Firstly, expand from the design data the ratio of clear area to one pixel; then, obtain 
expanded data. This expanded data is such that the ratio of clear area to one pixel is obtained 
in binary value or multi-level value. Employing multilevel value makes it possible to more 
precisely represent the expanded data required. 

[0028] At this time, in the case of a chromium mask or single-layered phase-shift mask, 
expand design data of a single layer. Letting the transmission of the glass substrate of a 
single-layer mask be 1, the intensity transmission of an absorbing member or "absorber" be t 2 , 
and the phase difference be 4>, the absorber's resulting complex amplitude transmission factor 
or "transmission" becomes equal to t-exp(icf)) with i=V-l. At this time, letting the expanded 
data of such a mask (ratio of clear area to a pixel) be 0<E(m,n)< 1, the complex amplitude 
transmission distribution F(m,n) is given as: 

F(m,n) = E(m,n) + t-exp(icj))(l-E(m,n)) 



[0029] Alternatively, for a multi-layer mask such as the ternary or tri-tone mask and the 
Levenson mask, expand from the design data of each layer the ratio of clear area to one pixel; 
then obtain expanded data. This expanded data is generated such that the ratio of clear area to 
one pixel is obtained in binary or multi-level value. Using the multi-level value enables more 
precise representation of the expanded data involved. 

[0030] In the case of a ternary or tri-tone mask, letting the transmission of glass substrate be 
set at 1, the intensity transmission of an absorber be t 2 , and a phase difference be (J) with 
i^Z-1, the absorber's complex amplitude transmission is equal to t-exp(icj)). At this time, 
letting the mask's expanded data be 0<E(m,n)< 1 while letting expanded data of a half-tone 
pattern be 0<E'(m,n)<l, the complex amplitude transmission distribution F(m,n) is 
represented by: 

F(m,n) = E(m,n) + t-exp(i<|))-E'(m,n) 
[0031] In the case of the Levenson mask, letting the transmission of the glass substrate be set 
at 1.0, the intensity transmission of an absorber be 1.0, and the phase difference of the shifter 
be 4> with i=\/-l, the absorber's complex amplitude transmission becomes t-exp(i4>). At this 
time, letting the mask's expanded data be 0<E(m,n)< 1 while letting the shifter pattern's 
expanded data be 0<E'(m,n)<l, the resultant complex amplitude transmission distribution 
F(m,n) is represented as: 

F(m,n) = E(m,n) + exp(ic|))- E'(m,n) 
[0032] From the foregoing, it can be appreciated that the complex amplitude transmission 
distribution data is obtained through generalization with respect to the above-stated three 
different types of masks by employing the following equation: 

F(m,n) = E(m,n) + t • exp(icj))- E'(m,n) 
[0033] Next, let the complex amplitude transmission distribution data thus obtained in the 
way stated supra be input to the reference data calculator unit 7. 
[0034] Processing for calculation of this reference data will be explained below. 
[0035] The image intensity distribution upon irradiation of light, I(x 1 ,y i ), is represented by 
the following Equation (1). Here, the term B 0 (x,y) is the phase coherency coefficient, K(x,y) 
is an image of point-source light, and F(x,y) is the complex amplitude transmission 
distribution of an object of interest. The phase coherency coefficient B 0 (x,y) is represented 
by Equation (2), whereas the image of point source light K(x,y) is by Equation (3). 
Additionally, R* is a conjugate complex number. 



I<X,yi) = J|B 0 (x 0 - x' 0 , y 0 - y' 0 )F(x 0 , y 0 )F*(x' o , y' 0 ) 

K ( x i - x o, Yi - y 0 )K(x, - x' 0 , y ; - y , o )dx o dy 0 dx' 0 dy' o (1) 



B 0 (x r y) = 2J a (u)/u f u = 2na^x 2 + y 2 (2) 

5 4) is the ratio, NAJNA 0 of the numeral aperture NA C of condenser lens to the numeral 

aperture NA 0 of objective lens, X is the wavelength of light source, and Jj(x) is the first 
Bessel function of the first kind. Note here that Cartesian coordinates are normalized so that 
x*= NA 0X /A, y 1 = NA oy A. 

K(xy) = exp(-i27i-^) f f exp[C7i{&c 1 +T]y' +0.5z'-(e+r| 2 ) +^,n)}]^rfr| ^ 

where, NA 0 is the numerical aperture of an objective lens, X is the wavelength of light source, 
1 0 wherein normalization is done so that x - NA 0X /A, y' = NA oy /A and z - NA 0X /A, respectively, z 
denotes the degree of defocus, and W(£, n) denotes the wave front aberration of the objective 
lens. 

[0036] Next, digitizing the Equation (1) on a pixel-wise basis, we obtain Equation (4) as: 
I(m,n) = £ZZX B 0 (mo - m' 0 , n 0 - n' 0 )F(m 0 , n 0 )F*(m , 0 , n' 0 ) 
1 5 K(m - m 0 , - n 0 )K*(m - m' 0 , n - n' 0 ) (4) 

[0037] Note here that although any one of K(m,n) and B 0 (m,n) has a maximum value at the 
origin, each of them rapidly becomes diminished with an increase in distance from the origin. 
Accordingly, the image intensity distribution I(m,n) may be approximated by Equation (5) 
given below: 

20 I(m, n) = B 0 (0,0) • { [ K(m, n) | 2 + | F(m, n) | 2 } 

+ 2 • Re[F*(m, n) • K*(0,0) • {(B 0 (m, n) • K(m, n))° * F(m, n))}] (5) 

[0038] The first term indicates the case of (m o -m' 0 , n o -n' o )=(0,0) whereas the second term 
shows the case of establishment of (m 0 -m' 0 , n o -n' o )*(0,0) and (m-m 0 , n-n o )=(0,0) or, 
alternatively, (m 0 -m' 0 , n o -n' o )^(0,0) and (m-m' 0 , n-n' o )=(0,0). Here, "*" is used to denote a 
25 convolution integral. To be brief, we obtain f(m,n)*g(m,n) = ££f(m-m 0 , n-n 0 )-g(m 0 , n 0 ). 
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[0039] In addition, Equation (4) is formulated into Equation (6) so that the image intensity 
distribution I(m,n) is represented by complex coefficients. At this time, P(m,n) is a finite 
impulse response (FIR) filter comprised of real coefficients whereas Q(m,n) is a FIR filter 
formed of complex coefficients. 
5 I(m, n) = P(m, n)* | F (m, n) | 2 + Re[F * (m, n){Q(m, n) * F(m, n)}] (6) 

where Re[X] represents the real part of a complex number X, and in the case where 
P(m,n) = B o (0,0) | -IK(m,n) | 2 , 
Q(m,n)= 0 in case of(m,n) = (0,0), 
and Q(m,n)= 2K*(0,0)-B0(m,n>K(m,n) otherwise. 
1 0 Additionally, when setting the complex coefficient FIR filter 

Q(m,n)-Q r (m,n)+iQ j (m,n), F(m,n)=F T (m,n)+iF i (m,n), we obtain from Equation (6) the 
following Equation (7): 

I(m, n) = P(m, n)* {F r (m, n) 2 + F ; (m, n)2} 

+ F r (m, n) {Q r (m, n)*F r (m, n) - Q;(m, n)* F^m, n)} (7) 
1 5- + F ; (m, n) {Q r (m, n)* F^m, n) + Q^m, n)* F r (m, n)} 

[0040] As apparent from Equation (7), by setting the complex coefficient of the FIR filter 
Q(m,n) in the way noted above, the image intensity distribution I(m,n) is indicated by 
addition of five real parts as will be shown below. 

I (m,n)=P(m,n)*{F r (m,n) 2 + Q r (m,n) 2 } + F r (m,n)Q r (m,n)*F r (m,n)-F r (m,n)Q i (m, n)*^ 
20 (m,n) +F i (m,n)Q r (m,n)*F i (m,n)-F i (m,n)Q i (m,n)*F r (m,n) (8) 

[0041] Thus it is possible by using the FIR filter to obtain from the design data the reference 
data that correspond to the scanned image data. 

[0042] An explanation will next be given of one practical implementation of the processing 
method for obtaining the complex coefficient FIR filter by the five real FIR calculation 
25 method discussed above to thereby calculate reference data (image intensity distribution 
I(m,n)), with reference to Fig. 3. 

[0043] The example shown herein is an example that calculates the reference data from 
design data of a specific mask structure with a mixture of a chromium mask and a phase-shift 
mask, known as a ternary or tri-tone mask among those skilled in the art. A first layer 
30 indicates clear patterns; a second layer shows phase-shift mask patterns. Note here that any 
overlap occurring between the first and second layer patterns on the design data should be 
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required to be excluded. Overlap is removed from test data by data processing prior to the 
inspection procedure. Alternatively, such overlap removal is also attainable by excluding an 
overlap on a real time basis through logical arithmetic processing from input data of the first 
and second layers, prior to obtaining the complex transmission distribution data and 
5 alternatively complex reflection distribution data. 

[0044] Letting the complex amplitude transmission of the phase-shift mask be t-exp(icj)) as 
stated previously, an appropriate value that corresponds to t-cos((|)) is set as a coefficient 1 at 
a coefficient hold/storage unit 13. In addition, a value corresponding to t-sin(4>) is preset as 
coefficient 2 at a coefficient hold/storage unit 14. 

10 [0045] As shown in Fig. 3, the first layer's expanded data is input to and stored in a register 
10. The second layer's expanded data is input for storage to a register 9. Next, the second 
layer's expanded data is input to a multiplier 100 for multiplication with the coefficient 1 as 
output from the coefficient storage unit 13 and then input to an adder 1 1 . Substantially 
;;i simultaneously, the first expanded data is also input to the adder 1 1 to add to the 

15 multiplication result of the second expanded data and coefficient 1. In this way the real part 
F T is thus obtained, which will then be input to a register 103. 

; ■ j 

|ll [0046] On the other hand, the second layer's expanded data as output to a multiplier 101 is 

~ multiplied with the coefficient 2 as output from the coefficient storage unit 14 to thereby 

• obtain imaginary part F i5 which will be input to a register 1 04. 
20 : [0047] Next, the real part F r as an output from the register 103 is input to a line buffer 20 and 

: - • then stored therein. The imaginary part F ; being an output from the register 1 04 is input to a 
line buffer 21 and held therein. The real part F r being held in the line buffer 20 and the 
imaginary part F; as held in the line buffer 21 are input simultaneously to an image intensity 
processing unit 15 in conformity with the significance of the FIR filter used. This image 
25 intensity processor unit 1 5 uses the real part F r and imaginary part F { to obtain a value of 
F r 2 +Fi 2 , which is then output to a data hold/storage unit 105. 

[0048] Additionally, the real part as output from the line buffer 20 is input to a register 106 
and held therein, whereas the imaginary part F 1 as output from the line buffer 21 is input to a 
register 107 for storage therein. 
30 [0049] Next, the F r 2 +F; 2 data as output from the register 105 is input to the FIR filter 16 for 
convolution integration with the FIR filter P of the real number coefficient to thereby obtain 
the result of P*(F r 2 +F 3 2 ) data (first term of Equation (7)), which is then output to register 108. 
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Note that the symbol "*" represents the convolution integral operation. 
[0050] Meanwhile, the real part F r , being held in the line buffer 20, is input to the register 
106 and stored therein. The imaginary part F ; , as held in the line buffer 21, is input for storage 
to the register 107. 

5 [0051] The real part F r , being held in the register 106, is input to a complex FIR filter 27. 

This complex FIR filter 27 effectuates convolution integration of Q r for calculation of F r *Q r 
and then outputs it to a subtracter 17. The imaginary part F ; , being held in the register 107, is 
input to a complex FIR filter 29. This complex FIR filter 29 performs convolution 
integration of Q, to calculate F*Q, and then outputs to the subtracter 17. The subtraction 
10 result of the subtracter 1 7, i.e., F^C^-F^Q, data, is then output to a multiplier 111. The real 
part F„ as output via a delay circuit 18, will also be output to the multiplier 111. This 
multiplication result F^F^Q.-F^Qj) data (second term of Equation (7)) is input to register 
109 and then retained therein. 

[0052] The real part F r , being held in the register 106, is input to a complex linear response 
15= filter 28. The complex FIR filter 28 convolution-integrates Q ; to calculate F r Q ; and then 
outputs to an adder 19. The imaginary part F i5 as held in the register 107, is input to a 
complex FIR filter 30. This complex FIR filter 30 applies convolution integration to Q T for 
calculation of and then outputs to the adder 19. The addition result of adder 19, i.e., 
F^Qj+F^Q,., is output to multiplier 112. The imaginary part F i? as output via the delay circuit 
20 1 8, is also output to the multiplier 1 12. Its multiplication result F/F^Q^F^Q,.) data (third 
term of Equation (7)) is input to a register 110 and then stored therein. 
[0053] In the way stated above, the result of the processing executed by the complex FIR 
filters 27-30 is such that the real part is given by Q^F^Q^Fi whereas the imaginary part is 
given as F^Qj+F^Q,. Real part output is achievable because of the fact that the real part 
25 Q r *F r -Q I *F i is multiplied by real part F r while letting the imaginary part F r *Q +Fj*Q r be 
multiplied by imaginary part F ; . 

[0054] Further, the data items (second and third terms of Equation (7)), being held in the 
registers 109, 110, are added together at an adder 26, causing its addition result to be added at 
an adder 25 to the data (first term of Equation (7)), as held in the register 108, to thereby 
30 obtain a final addition result for use as an image intensity distribution (right part of Equation 
(7)), which is then output as the reference data required. 

[0055] The reference data thus obtained is input to the comparator unit 8 shown in Fig. 2 and 
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will then be subject to comparison with scanned image data. 

[0056] Turning to Figs. 4A to 4D, there are shown some major image patterns each 
indicating the image intensity of the reference data as obtained in the method stated supra. 
Figs. 4A and 4C show optical images of the bright field whereas Figs. 4B and 4D show 
optical images of dark field. Additionally Figs. 4A-4B each depict the reference data as 
obtained through filter calculation in accordance with this embodiment. Figs. 4C-4D are 
simulation results, each of which was calculated while applying a rigorous partial coherent 
model to an identical pattern. 

[0057] It is readily apparent, by comparison of the bright field cases of Figs. 4A and 4C, that 
the pattern of Fig. 4A, obtained through calculation in consideration of the phase of light 
using the FIR filter of complex coefficients, is identical with maximal fidelity to the pattern 
of Fig. 4C, which was calculated using the rigorous partial coherent model. These figures 
demonstrate that the use of such a scheme for obtaining optical image(s) by calculation in 
view of light phase, using more than one complex coefficient FIR filters that are efficient in 
this calculation method, may result in successful representation of the difference in image 
profile between bright and dark fields, and also the difference in image profile between the 
chromium mask and phase-shift mask, which differences have been difficult to represent in 
the prior art. 

[0058] Similarly, with the dark field cases of Figs. 4B and 4D also, it is apparent that the 
pattern of Fig. 4B, as calculated in view of light phase using the FIR filter of complex 
coefficients, is nearly identical to that of Fig. 4D as calculated by the strict partial coherent 
model. 

[0059] In this way, in accordance with the present invention, it becomes possible to obtain 
by simulation the required reference data, which is nearly identical to the optical image of a 
test substrate as has been sensed or picked up upon actual irradiation of light thereto, which 
in turn makes it possible to preclude creation of any unwanted recognition errors of the 
inspection apparatus used, thus enabling further improvement in accuracy of die-to-database 
schemes. 

[0060] It has been stated that, according to the present invention, high-accuracy optical 
images may be successfully simulated from design data, even in those masks employing 
resolution enhancement techniques; thus it is possible to provide an apparatus which prevents 
nuisance defects from being detected, even when using it as reference data for comparison. It 
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is also possible to calculate the reference data with maximal accuracy in accordance with 
optical illumination conditions (light source wavelength of inspection apparatus, numerical 
aperture (NA) of an objective lens, partial coherence factor, degree of defocus or the like), 
lens aberration (spherical, astigmatism, coma aberrations or otherwise), mask absorber's 
physical properties (transmission amplitude, reflection amplitude, phase differences at the 
inspection wavelength), and others. 

[0061] In the description above, the reference data was obtained with respect to all surfaces 
of the test substrate in an all-at-a-time fashion for comparison with the optical image of a 
specimen or sample. The present invention should not be limited only to this method, and 
may also be arranged so that a pixel immediately before the to-be-compared pixel is obtained 
through simulation as the reference data and then subjected to comparison with an optical 
image corresponding to such pixel. 

[0062] Additional advantages and modifications will readily occur to those skilled in the art. 
Therefore, the invention in its broader aspects is not limited to the specific details and 
representative embodiments shown and described herein. Accordingly, various modifications 
may be made without departing from the spirit or scope of the general inventive concept as 
defined by the appended claims and their equivalents. 
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